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ABSTRACT 

Mars is likely to be a planetary embryo formed through collisions with plan- 
etesimals, which can explain its small mass and rapid formation timescale ob- 
tained from ^^^Hf-^^^W chronometry. In the classical theory of planet formation, 
the final embryo mass is determined only by the solid surface density. However, 
embryos can stir surrounding planetesimals, leading to fragmentation through 
erosive (cratering) collisions. We find that radial drift of small fragments can 
drastically reduce the solid surface density. On the other hand, embryo growth 
is accelerated by fragment accretion. Since collisional fragmentation efficiency 
depends on the initial size of planetesimals, the final embryo mass and its growth 
time are determined by the initial planetesimal size and disk surface density. We 
have investigated the effect of these two parameters on the mass of Mars and the 
predicted radiogenic excess of ^^^W in the martian mantle. Two scenarios can 
explain the rapid formation of small Mars: (i) it formed by accretion of small 
planetesimals in a massive disk or (ii) it formed from large planetesimals but its 
growth was arrested by the inward then outward migration of Jupiter. Taking 
into account all constraints, we conclude that Mars is likely to have formed in a 
massive disk of about ~ 0.1 solar mass from planetesimals smaller than ~ 10 km 
in radius. Such small planetesimal size cannot explain core accretion of Jupiter, 
suggesting that there may have been a heliocentric gradient in planetesimal size 
in the solar nebula. 

Key Words: Mars; terrestrial planets; planetary formation; planetesimals; origin, solar 
system 
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INTRODUCTION 



Mars, one of four terrestrial planets in the solar system, has a small mass of 0.11 
Mq (M^ =Earth mass) and an orbital semimajor axis of 1.5 AU. Terrestrial planets 
are believed to have formed via collisions among Mars-sized planetary embryos in the 



long-term orbita 



Chambers et al. 



2006 



instability that fo 



1996; 



Ogihara et al. 



Iwasaki et al. 



owed substantial gas dissipation (i.e., chaotic growth ; 



2001 



Kominami and Ida 



2002; 



Kenyon and Bromley 



20071 1. The ^^^Ri-^^ ^W short-hved chronometer (t 



very short accretion timescale for Mars (JDauphas and Pourmand 



dissipation time in 



2011 



Haisch et al 



erred from disk and meteorite observations (JDauphas and Chaussidon 



/2 = 9 Myr) gives a 



20111). similar to the gas 



200 ll ). This indicates that Mars is likely to be a remnant of the early 



generation of planetary embryos formed via collisional accretion of planetesimals. 

Collisions between planetesimals produced planetar y embryos, which grew fur- 



ther through collisions with surrounding planetes imals (iWetherill and Stewart 

3 



1993 



Weidenschilling et al. 



1991 



Kokubo and Ida 



19981). Once embryos became massive, they 



stirred surround ing planetesimals, producing collisional fragme i its that contributed to 



embryo growth fiWet 



lerill and Stewart 



Kenyon and Bromley 



2004 



2OO9I : 



1993 



Inaba et a. 



Ormel and Kobayashi 



2003 



Chambers! i2006 



2008 : 



2012). When the surrounding 



bodies were depleted due to the radial drift of small fragment s and/or accretion onto 



embryos, the growth of embryos stalled (jKobayashi et al. 



2OIOI ). In oligarchic growth, the 



initial size of planetesimals and disk mass influence the final masses of embryos and their 
accretion timescales. Small planetesimals are easily fragmented during collisions. The 
resulting fragments are small enough to feel some drag from the gas, leading them to spiral 
towards the Sun. As a result, small planetesimals tend to produce embryos that have 
smaller masses than when large initial planetesimals are considered. On the other hand, 
the interaction of fragments with the gas damps their relative velocities and increases the 
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efficiency with which they can be accreted by embryos through gravitational focusing. 
Thus, with small planetesimals the timescale of embryo growth is shorter than when large 
planetesimals are considered. 

During formation of the terrestrial planets, the segregation of metal from silicate 
produced metallic cores, the timing of which can be estimated using ^^^Hf-^^^W systematics. 
Hafnium, a lithophile (i.e., rock-loving) element, was retained. In contrast, tungsten, a 
siderophile (i.e., metal-loving) element, was partitioned into metal cores. Hafnium-182 
(i^^Hf) decays to tungsten-182 {^^"^W) with a half life of 9Myr. The decay of ^^^^Hf in 
the martian mantle after removal of W in the core resulted in the production of excess 
radiogenic ^^^W compared to other isotopes of W in chondrites. The ^^^W excess in the 
martian mantle is known from measurements of shergottite-nakhlite-chassigni te (SNC) 
mete o rites t hat are thought to have been ejected from Mars b y impactors {e.g 



2004 



2009 



Foley et al. 



2005 : 



Dauphas and Pourmand 



Kleine et al. 



20111). Strictly speaking, ^^^W 



variations record the timing of core formation but this process presumably tracked planetary 
accretion. In this contribution, we use ^*^W excess in the martian mantle and the mass of 
Mars to constrain the disk surface density and initial planetesimal size in the Mars-forming 
region. 

2. MARS FORMATION 



We use the following disk model for the initial surface mass density of planetesimals 

Ss,o and gas Sg^o^ 

a \-3/2 



s,0 



-'gas,0 



a;SMMSN,g 



gem ^, 



lAU. 
a \-3/2 

lAuJ S'"^" ' 



(1) 
(2) 
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where S 



MMSN,s 



7.1 Kcm ^ and S 



MMSN,g 



1.7 X 10^ gem ^ are, respeetively, the solid and 



gas s urface densities at 1 AU in the minimum- mass solar nebula (MMSN) model (JHayashi 



198ll ). and a is the distance from the Sun. We vary the scaling factor x to investigate the 
conditions for Mars formation. The solid surface density evolves naturally in the simulation, 
while the gas density is artificially imposed to decrease as a function of time as exp(— t/ra) 
with a gas dissipatio n timescale ta. Dis k observations indicate that r^ is on the order of 



several million years (JHaisch et al. 



20011). 



Kobayashi et al.l (120101 l201l[ ) showed that model outputs of statistical simulations 
of runaway and oligarchic growth are controlled primarily by the planetesimals that 
dominate the initial solid surface density, so the results are largely insensitive to the 
initial size distribution of other minor bodies. For this reason, we have adopted a fixed 
initial planetesimal radius rn rather than using a n actual initial size d istribution, which 



is poorly known anyhow ( iMorbidelli et al 



2009 



Weidenschilling 



20111 ). We limit our 



analysis to tq > 1km because in a turbulent disk, smaller planetesimals would not be 



sensitive to gravitational 



Youdin and Lithwick 



2007; 



bcusing, thus preventing runaway growth to proceed (e.g. 



Ormel and Cuzzi 



2007; 



Ormel et al. 



20101). 



We performed numerical simulations of Mars formation via collisional evolution of 
planetesimals. In the model, the disk is divided into four annuli. The width of each annulus 
is given by 0.2 times the annulus' characteristic radius defined by the mean of the inner and 
outer r adii of the annulu s . The characteristic radii of the annuli are set to 1.5, 1.8, 2.2, and 
2.7AU. 



Kobayashi et al 



(120101, 



201 if ) showed that because the timescale for dust removal 
by gas drag is much shorter than the timescale for embryo growth, considering a larger 
disk did not change significantly the final embryo mass. The mass distribution of bodies in 
each annulus is treated as discrete mass batches with the mass ratio between adjacent mass 
batches fixed at 1.2. We follow the evolution of the mass distribution via collisions between 
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bodies. The coUisional rates depend on the relative velocities between colliders, which 
depend on their eccentricities and inclinations. The mass distribution affects eccentricities 
and inclinations by collisionless interaction between bodies. We the refore calculate th e mass 
and velocity evolution simultaneously, applying the collision rate s (llnaba et al.l 120011 ) and 



the velocity evol ution rates due to mutual interactions of bodies ( jOhtsuki et al. 



20021) and 



due to gas drag (lAdachi et al.iil976i ). We do not treat collisions and collisionless interactions 
between bodies in different annuli and therefore adopt an annulus width larger than the 
feeding zones of planetary embryos. The total fragment mass ejected by a single impact 

rn-i + m'))(t>/(l + (t>) using th e 



between colliders with masses mi and m2 is assumed to be 

dimensionless impact energy = mim2V^/2{mi + m2)^QD (jKobayashi and Tanaka 



2010 



Kobayashi et al. 



20101 ). where v is the coUisional velocity between the colliders and Qd is the 
specific impact energy (ergg~^ of impacted material) needed for the ejection of a half-mass 
of colliders. Hydrodynarn ic impact simulations give QJj, for 100 km sized or smaller bodies 



( JBenz and Asphaug 



19991 ). while 



the gravitational binding energy (jStewart and Leinhardt 



Qn of bodies with radius r > 100 km are determined by 



20091 ) . Connecting the two regimes 



(above and below about 100 km), we express Q^ value of a body with mass m, radius 
r, and density p as Q^ = (5os(^/l cm)^" + Qogp^r /Icm)'^'^ + 2CggGm/r, where G is the 
gravitational constant and C^g = 9 from coUisional simulations with distinct element method 

20091 ). Here we use Qos = 3.5 x lO^ergg"^, Qog = 0.3ergcm^g~^, 



(iStewart and Leinhardt 



Ps = —0.38, and (3 ^ = 1.36, which were obta ined from hv drodynamic simulations 



and 



Inaba et al. 



wit h V = 3km/s j Benz and Asphaug 



19991). Following 



Wetherill and Stewart 



or basalt 



fll993h 



(J200lh . when the bodies reach a certain mass mj-^n su ch that the sum o f 



their mutual Hill radii equals the radial-bin width divided by 6 = 10 ( jKokubo and Ida 
20021 ). the bodies are regarded as "runaway bodies" that do not experience collisions and 
dynamical interactions due to close encounters among them. Using statistical methods. 



we can accurately calculate embryo formation and growth until the subsequent stage of 
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chaotic growth, which may not be significant for Mars. The statistical simulation used 
in the present contributio n reproduces the resu 



fragmentation is ignored (JKobayashi et al. 



ts of A^-body simulations when coUisional 
20101) ■ It also reproduces the analytical so lution 



for mass loss due to a collisional fragmentation cascade (JKobayashi and Tanaka 



2010h . 



Embryos grow via collisions with surrounding planetesimals and fragments. If we 
assume that the bodies surrounding embryos are not lost at all, then embryos can accrete all 
the bodies in their feeding zones. The maximum mass of the embryo is called the isolation 



mass and is given by (e.g.. iKokubo and Idall2002l ) 



M,o = 0.13x3/^ 



^« 



(3) 



.1.5AU. 
where we use Eq. ([1]) for the solid surface density. Therefore, a disk with x ~ 1 produces 

Mars-size embryos at 1.5 AU if embryos can collect all solids. However, since small 

planetesimals are vulnerable due to low gravity, mass loss can occur by collisional 

fragmentation of planetesimals with radial drift of the resulting fragments, so the final 

embryo mass can be smaller than Miso. Destructive collisions between planetesimals 

produce small bodies, which can in turn be broken apart resulting in still smaller bodies. In 

the collisional fragmentation cascade, the mass distribution of bodies follows a power law; 

the surface number density of bodies with mass ranging fro m m to m + dm is proportiona l 



2010l ). 



to m~°'dm, where a = (11 + 3p)/{Q + 3p) with v'^/Q^ oc m^ ( JKobayashi and Tanakal 
The mass distribution is determined by the collisional cascade for bodies with r < ro, while 
runaway growth of planetesimals regulates the mass distribution of bodies with r > tq. 
With such a mass distribution, planetesimals with radii ~ ro determine the solid surface 
density. The timescale of planetesimal depletion by the collisional fra gmentation depends 



on Q^ of planetesimals and is given by (JKobayashi and Tanaka 



2010|) 



l.lxl0^x-°-^if^^ 
V 10 km 

-0.48 



1.69 



a 



1.5 AU 



3.2 



X 



M 



0.11 Mff 



Qog 



O.Sergcm'^g 



-0.72 

3.-2/ y^' 



(4) 



where we assume that the planetesiinal eccentricities are in an equihbrium state controlled 
by stirring by embryos and damping by gas drag and use Qjj = QogPl'^o/l cm)''^ vi^ith 
p = 3gcm~'^ for planetesimals. Erosive (cratering) collisions with spe cific impact energies 



smaller than Q^ control this timescale (JKobayashi and Tanaka 



20101 ). which is 4-5 times 



shorter than that derived b y ignoring erosive collisions as previous studies did (e.g., 



Wetherill and Stewart 



19931 ). This timescale was obtained by assuming power law mass 



distribution in the collision cascade. Several effects can lead to departure from such 
power-law behavior but this would not change significantly the timescale because it is 
determined by bodies for which the mass distribution does not deviate considerably from 
a power-law dist ribution with a ~ 11/6. Sin ce t^c is inversely proportional to the solid 



surface density (JKobayashi and Tanaka 



2010l ). Tec increases with decreasing Eg and the 



full-depletion timescale of planetesimals is much longer than Tc, 



The timesca^ 



2002 



Chambers 



e for embryo growth via plan etesimal accretion is given by (JKokubo and Ida 



2006 



Kobayashi et al 



2010), 



Tr 



grow,p 



2.7 X 10^ x"^/^ (- 



X 



ro_\ 
lOkmy 
a \ 27/10 / M 

l.SAuJ 



2/5 



1/3 



yr. 



(5) 



' grOW,p r^ ''"c(; . 



^0.11 M( 
For To > 100 km, embryos can reach Miso in a disk with x ~ 1 because r^ 
However, the final embryo mass is smaller than the isolation mass for tq < 100 km. Larger 
X is needed for smaller planetesimals to compensate for loss of mass through gas drift of 
collision fragments and achieve the same final embryo mass. Note that the necessary x to 
produce Mars-size embryos is lower than x given by pos ing Trr = r^mw^p, becau se Tcc and 



'^grow,p increase when the solid surface density decreases (JKobayashi et al. 



20101). 



For small values of vq (< 1-10 km), embryo growth proceeds to a large extent by 
fragment accretion. The coUisional fragmentation cascade halts when eccentricities and 
inclinations of the bodies are damped well by gas drag. Fragments thus accumulate at 



-9- 



' ~ 0.5Qd. This condition gives a typical fragment size of 1-10 m depending on the 
embryo mass M. The radial drift timescale of such bodies is given by (see the derivation in 



Kobayashi et al. 



2OIOI) 



r. 



'''drift 



stop 






6.9 X 10' 



1.5 AU 

Qh 



M 



6.1 X lO^ergg^i 



0.11 Mg 



(6) 



where Qk is the Keplarian frequency, we use the value of Q^ for r = 1 m, Xstop is the 
stopping time by gas drag for the fragments under the Stokes regime, vk is the Keplarian 
velocity, Pvs,iow = 73 is the dimensionless viscous stirring rate for low eccentricity, and 
V = {vk — v^as)/vK- Since small fragments have low eccentricities and inclinations, fragment 
accretion allows fast growth of embryos. T he embryo growth timescale via fragment 



accretion is given by ( JKobayashi et al 



201 



'grow, frag 



4.3 X 10^ x-^ 



M 



0.11 Mff 



1/3 



1.5 AU 



yr. 



(7) 



If we ignore the radial drift of fragments, e mbryos reach Miso in a short timescale of Tgrow.frag 



(JRafikov 



200 



Kenyon and Bromley 



20091 ). However, fragments are removed by radial drift 



due to gas drag in a shorter timescale of Thrift, which stalls embryo growth before they reach 
isolation masses. In the case of fragment accretion (small initial planetesimals) , formation 
of Mars-size embryo requires Tdrift ~ Tgrow.frag; 2; > 6. With larger initial planetesimals, 
coUisional fragmentation of planetesimals is less effective and therefore embryos can reach 
the mass of Mars for smaller x. 

To summarize, small initial planetesimals are prone to fragmentation and mass loss 
occurs through gas-induced drift of fragments. As a result, small planetesimals tend to 
produce small embryos (and vice versa). The eccentricities and inclinations of fragments 
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produced by collisions between small planetesimals are damped by interaction with the 
gas. The resulting reduction in relative velocities increases the gravitational cross section of 
embryos for fragments, so small planetesimals allow rapid embryo growth (and vice versa). 

We start the investigation for embryo growth in a MMSN disk, where the final embryo 
mass is about the mass of Mars if embryos can fully accrete bodies in their feeding zones; 
Miso ~ O.IM® at 1.5 AU. Fig. [1^ shows the time evolution of planetary embryo masses at 
1.5 AU for ro = 1, 3.7, 27 and 100 km with t^ = lOMyr. Growth rates and final masses 
of embryos depend on tq. Mars-sized embryos can be formed from planetesimals with 
ro > 10 km in the MMSN disk, while a more massive disk is necessary to form Mars from 
smaller initial planetesimals. 

The embryo mass evolution affects the predicted ^^^W radiogenic excess caused by 
the partitioning of W into the core and the decay of ^^^Hf in the mantle. In the embryo 
mantle, the excess abundance of ^^^W compared to the chondritic reference val ue (CHUR 



ch ondritic uniform reservoir) is g iven by (see the derivation in Appendix lAj 



2005 



Dauphas and Pourmand 



Jacobsen 



201 ih 



pl82.,r _ ^ / ^ .Hf/W , 

fc VV mantle — VW [ igOJJf I ^mantle ^ 



l82Hf^^ 

J mantle 

Hf/W 



X f (^^) " exp(-At')rft', 



/o \M{t))J 

where e'^'W^.,^,,, = [e''W/'''WU^,iJC''W/'''W)cnvR - 1] x 10^ gw = 
(i80Hf/i82w)cHUR X 10^ f^'JZ = (Hf/W)mantie/(Hf/W)cHUR - 1, M(t) is the em- 
bryo mass at time t, and A is the decay constant of ^^^Hf. From previous studies, 
^'^'W^Marsmantie = 2.68 ± 0.19, f S ^ ntb = 3-38 ± 0-56, Qw = 1.07 X 10^ and 



n82Tjf /180 



Hf/i8°Hf) = (9.72 ±0.44) x 10"^ fiDauphas and Pourmand 



2011 



and references therein) . 



Our simulations generate as model output the time evolution of an embryo mass M(t), 
allowing us to calculate s^^'^Wraa.ntie- Fig. [T]d shows the evolution of e^^^Wmantie during 



- 11 - 



embryo growth with different initial planetesimal radii in the MMSN disk. The predicted 
e^^^WMars mantle values are smaller than the measured one for r^ = 1-100 km. Small tq values 
produce e^^^Wmantie close to the martian mantle value (Fig. [Dd) but the final embryo mass 
is too small and a more massive disk is required (Fig. [1^). 

The results for massive disks with x = 2.7 are shown in Fig. HJ Smaller initial 
planetesimals with ro = 3.7km can make Mars-size embryos, while larger planetesimals 
make embryos that are much larger than Mars. The predicted e^^^W value for ro = 3.7km 
is about 2.84, which matches the observation. Therefore, a parameter set of x ~ 3 and 
ro ~ 4km can account for the mass and e^^^W value of Mars. For x = 2.7 and ro = 3.7km, 
the two inner annuli produce 8 Mars mass bodies, while smaller bodies are formed in the 
outer annuli close to the outer edge of the simulation. The predicted e^^^W values of the 
embryos are similar but slightly lower for embryos formed in the outer disk relative to those 
formed in the inner disk. In the proposed scenario, one of the embryos becomes Mars while 
the others accumulate to make Earth and Venus size planets in the subsequent stage of 
chaotic growth (see §3.2p . 

To investigate the conditions required for Mars formation, we have evaluated 
systematically the effects of initial planetesimal size and solar disk mass on the final mass 
and accretion timescale of embryos in the Mars-forming region. Fig. [3^ summarizes the 
final embryo masses and the ^^^^VFinantie values derived from 88 simulations at 1.5 AU for a 
logarithmic grid of a; x ro, where x = 1, 1.4, 1.9, 2.7, 3.7, 5.2, 7.2, 10 and ro = 1, 1.9, 3.7, 
7.2, 14, 27, 52, 100, 190, 370, 720 km. Final embryo masses reach the mass of Mars on the 
thick solid curve in the x and ro space of Fig. |3^. Our statistical simulations present a mass 
distribution of embryos that spans a factor 4-5 in mass, which is a measure of uncertainty in 
the predicted embryo mass (thin solid curves). The tendency that small initial planetesimals 
require a massive disk for the formation of Mars-sized bodies is explained by the analytical 



-12- 



solutions for the final embryo mass derived by iKobayashi et al.l ( 120101 ) . Mass loss due to 
collisional fragmentation with radial drift of fragments reduces the final mass for small tq. 
The e^^^Wmantie ± ^^^^^Wmantie values are calculated from Eq. (jHD, where Ae^'^^Wmantie is 
the error in the predicted e^^^Wmantie value calculated by propagating all uncertainties on 
model parameters. At 1.5 AU, the parameter space that reproduces both the accretion 
timescale and final mass of Mars is tq < 10 km and x > 1. 

To evaluate the robustness of our conclusion, we have carefully investigated the 
parameters that could affect the mass and accretion timescale of Mars. In particular, the 
influences of embryo atmospheres and gas drag on accretion dynamics, embryo migration, 
disk dissipation timescales, and disruption strength of bodies are discussed below. 

Fig. H^ shows the results of simulations including the collisional enhancement for 
capture of small particles by embryos induced by the presenc e of a planetary atmosp here, 



20031 ). The 



with a density profile determined by the atmospheric opacity (llnaba and Ikomal 
parameter space that explains the mass and e^^^W of Mars is similar to the case without 
the collisional enhancement by planetary atmosphere. The reason is that embryos have 
small atmospheres and the atmospheric effect is negligible for an embryo smaller than Mars 



teobayashi et al. 



20111 ). 



As discussed previously, gas can d amp velocities of small p articles, thus increasing the 
efficiency of their capture by embryos. lOrmel and Klahrl (J2010l ) further found that gas can 
directly affect the collisional cross-sections between embryos and small bodies. While the 
dissipative nature of gas drag promotes collision, strong particle-gas coupling acts against 
it by allowing fragments to be carried along with the gas passed the embryo. Fig. |3)d 



illustrates the results o: 
drag that is derived by 



simulations including t he collisional cross section affected by gas 



Ormel and Klahrl (l2010l ). The collisional enhancement is significant 



only at Tstop ~ ^k ^^^ Mars-size embryos. Since the typical fragments have radii ~ 10 m 
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corresponding to Tstop ^ ^k^ ( Kobayashi et al.ll2010l ). this coUisional enhancement is minor 
for Mars formation (compare Figs. [3^ withHb). 

We additionally investigated how the gas dissipation time Td affected the accretion 
history of Mars. The final embryo mass and e^^^W value for r^ = 5 Myr are almost the same 
as those for tj = 10 Myr (Fig. \51^). Fig. \5jp shows the results for an extremely short gas 
dissipation time of r^ = 1 Myr. The parameter space where the mass and accretion timescale 
of Mars are reproduced is not very different from that obtained when the dissipation time 
is set to higher values of 5 or 10 Myr. 

Recent studies suggested that terrestrial planets were c onstructed from a narrow ring 



betw' e en the current orb its of Earth and Venus (0.7 to 1 AU; iMorishima et al. 



2009 



Walsh et al. 



2011 



2008 



Hansen 



. In this scenario, Mars was likely scattered from around the outer 
edge of the truncated disk. We thus additionally investigate the possibility of Mars 
formation at 1 AU by setting the characteristic radii of the innermost and outermost annuli 
at 1 and 2.7AU, respectively. Figure [3b shows the final embryo mass and the e^^^W 
value at 1 AU in the tq and x space. The parameter space required for the production of 
Mars-sized embryo is similar to the case at 1.5 AU. However, larger planetesimals with 
tq < 40 km could explain the mass and accretion timescale of Mars. 



A very important 



2010 



develo pment of this work and other recent studies (JKobayashi and Tanaka 



Kobayashi et al 



20101 ) is the recognition that erosive collisions and fragment accretion 
dominate embryo gro wth. The treatrnent of fragmentation depends on Qp, w hich affects 



final embryo masses (JKobayashi et al 



2010 



2011 



Bromley and Kenyon 



201l|). The Q*^ 



value is uncertain. For instance, bodies m ade of rubble piles hav e Q^ smaller than what 



we used, thus producing smaller embryos (JKobayashi et al. 



20101 ). With ro < 10km, low 



Qd requires larger x for Mars formation, which may rise the value of plausible initial 
planetesimal radii only slightly. 
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3. Discussion 



3.1. Comparison with previous studies 



Terres t rial planet formation w as investigated by A^-body simulations (JKokubo and Ida 



1998 



2000 



Morishima et al 



2012I ). Our statistical simulation does not follow the orbits of 



bodies individually, but reproduces well the mass and random- veloc ity evolution obtainec 



from A^-body simulations under the assumption of perfect accretion (JKobayashi et al. 



2010|). 



For terrestrial planet formation, 



Wetherill and Stewart 



(119931 ) performed statistical 
simulations similar to ours and obtained embryos of mass ~ 3 x 10^^ g after about 
0.1 Myr for a; ~ 2 and tq ^ 7km, which is slightly larger th an embryo masses from our 



simulations (~ 2 x 10 ^ " g) 



Jsing the fragmentation model of iWetherill and Stewart 



(119931), 



Kenyon and Bromleyl ( 120041 ) carried out statistical simulations for vq < 1 km and obtained 
embryos of mass about 5 x 10^^ g at 1 Myr for x ~ 2. In our simulations, the embryo masses 
are of about 2 x 10^^ g at 1 Myr for x ~ 2 and tq = 1 km. Thus, the final embryo masses 
predicted in our simulations are systematically lower than those obtained by previous 
workers for the same sets of conditions. The difference stem s from a different treat r nent o f 
erosive (cratering) collisions. In the fragmentation model of IWetherill and Steward (19931), 



erosive collisions with specific impact energies smaller than Q^ yield less ejecta than those 
that are expected from impact laborat ory experiments, while these a re well reproduced by 



the fragmentation model developed by 
present contribution. 



Kobayashi and Tanaka 



(120101 ). which is used in the 



Weidenschillingl (J2008[ ) carried out simulations of terrestrial planet formation where 
planetary embryos were treated as individual entities, while small bodies were followed 
using mass bins. These simulations showed that smaller initial pl anetesimals produced 



smaller planetary embryos. This tendency is consistent with ours. 



ChambersI (]2008|) used 
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a similar mo d el; sta rting from planetary embryos and planetesimals. The simulations 
in IChambersI (|2008[ ) followed collisional fragmentation cascades using mass bins with a 



mass ratio between adjacent bins of 10. These wide mass bins artificiall y slowed down 



planet esim al depletion by th e fragmentation cascade by a factor of ~ 4. 



(120081) and 



Weidenschilling 



ChambersI (120081) considere d the migration of e mbryos due to tidal interaction 



with the nebula (type I migration; e.g., 



Tanaka et al. 



20021 ). We briefiy discuss the effect of 



type I migration on terrestrial planet growth below. 



3.2. Subsequent Evolution of Mars 



From our simulations, tens of Mars size embryos with orbital separation of 10 mutual 



Hill radii are expected to form in t he terrestrial planet for ming 



remain small during chaotic growth, iMorishima et al 



(120081) and 



region. 



Hansen 



^or Mars to 



fl2009[ ) suggested 



that planetary embryos must have been concentrated inside about 1 AU. Such orbital 
concentration may be attained by type I mig ration of embryos c aused by gravitational 



interaction with the surrounding nebula (e.g.. 



Tanaka et al. 



20021 ) if the torque cancellation 



due to the unsaturation of corotation torque occurs at < 1 AU w 



accre tion determines disk temperature (jPaardekooper et al. 



2010 



lere heating by disk 



2011 : 



Kretke and Lin 



20121 ). Mars may have evaded collisions with other embryos during chaotic growth by being 
scattered to its current orbit at 1.5 AU. The planets resulting from simul ations of chaotic 



grow th have orbital eccentricities higher than those of Venus and Earth (jChambers et al. 



19961 ). but these eccentricities were most likely da mped by gravitational in t eractions with 



a tenuous gas disk and/or a planetesimal swarm ( iKominami and Ida 



2006 



Morishima et al. 



2002 



O'Brien et al. 



20081 ). This dynamical friction could have been achieved by small 



< 10 m) leftov er planetesimals that w ere accreted by Earth and Mars after completion of 



core formation (jSchlichting et al 



20121 ). The mass of the late veneer on Mars is estimated 
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from abundances of highly siderophile elements in SNC me teorites to be ~ 1.6 x 10 g. 



Assuming a chondritic impactor composition 



this would have delivered ~ 1.6 x 10 g of W (iWarren 



100 ng 



-1. 



1999 



Wasson and Kallemeyn 



Walker 



Lodders and Fegley 



19881), 



2009 ). The martian 



19971 ). Thus, 



mantle contains ~ 3.2 x 10^^ g of W (62 ppb W; 
planetesimals responsible for damping the eccentricities and inclinations of the terrestrial 
planets may have delivered only ~ 0.5% of the W inventory of the martian mantle and 
would have had a negligible effect on the W isotopic composition. 



Walsh et al. 



( 120111 ) proposed that Jupiter migrated inward until Saturn was caught 
in a 2:3 resonance with Jupiter, at which point Jupiter migrated outwards. The U-turn 
migration of Jupiter at 1.5 AU would have truncated the disk and could explain Mars 
present orbit and mass. If the Jupiter's passing fully cleaned up small bodies, proto-Mars 
could no longer grow after it was scattered to the current orbit. Our simulations can be 
interpreted in light of this scenario by considering Mars' isolation time Tj, after which the 



201l|) 



embryo stops growing. The U-turn migration of Jupiter occurred at t^ (iWalsh et al. 
Proto-Mars could be scattered to its current orbit after the migration of Jupiter; n > t^. 
For Tj > 20Myr - several times longer than the gas dissipation time t^, Mars' small mass 
and rapid accretion is best explained by starting with small planetesimals as was discussed 
above (Fig. [6^). However, if one considers a rapid isolation t ime, which is more likely 



in the context of the Grand Tack scenario (jWalsh et al 



201 ll ). models that start with 



large planetesimals are allowed (Fig. |6]d), because Mars avoids late accretion of chondritic 
material. As discussed in Sect. 3.3, considering the conditions needed to form Jupiter and 
Saturn, the scenario of formation of Mars from large planetesimals in a MMSN disk is 
unlikely. 

A standing question is whether collisions between Mars and other embryos during 
chaotic growth could have lead to further growth and modified the W isotope evolution 
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of the martian mantle. During such colhsions, the core of the impactor may partially 
merge with the target core without complete isotopic equilibration with the target mantle 



( JDahl and Stevenson 



2010) ■ Therefore Mars growth via collisions of sub-Mars embryos 
in the subsequent stage of chaotic growth might be allowed to some extent. However, 
a planetesimal disk leading to small embryos would not have contained enough mass to 
produce the terrestrial planets, which makes it unlikely that Mars grew by the protracted 
accretion of small embryos. 

Mars has experienced at least one collision at a later stage. The northern hemisphere 
of Mars is lower and has thinner crust that the southern hemisphere, a feature that is 
commonly referred to as the martian hemispheric dichotomy. Several interpretations 
have been proposed, including one that calls for a large impact that would have 
obliterated the northern hemisphere less than 100 Myr after Solar System formation (e.g., 



Wilhelms and Sauvres 



20081). 



Marinova et al 



1984 



Andrews-Hanna et al 



2008; 



Marinova et al 



2008 



Nimmo et al. 



( 120081 ) used smooth particle hydrodynamics simulations to constrain 



the impactor to a 1,000 km-radius body that would have impacted Mars at a ~ 45° angle. 



Sotin et al. 



Using the mass-radius relationship R/R^ = (M/M^)^-^^^ (i?® is Earth radius; 
20071 ). such an impactor would have delivered ~ 2% of the mass of Mars. Most likely, 
this would have taken place after most ^®^Hf had decayed, possibly ~ 50 Myr after solar 
system formation. We have evaluated the impact of such a late impact on the W isotopic 
evolution of the martian mantle for one of the simulations that reproduces the mass and 
e^^^W values of Mars (x = 2.7 and tq = 3.7km). The calculated 5^^^W shifts from 2.7 to 
2.5 when taking into account this late impactor. This is well within uncertainties, so our 
conclusions are not affected by the possibility that the martian dichotomy wa s produced by 



a late impactor. The assumption of a 1000-km radius body is conservative as 



Nimmo et al. 



( 120081 ) considered a head on collision with a smaller body of 320 km radius, which would 
have a negligible effect on e^^^W. 
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3.3. Radial Gradient of planetesimal sizes 



As discussed in §3.2[ two scenarios can explain the small mass and e^^^W value of 
Mars: (i) small planetesimals in a massive disk formed Mars or (ii) large planetesimals 
in an ~MMSN disk formed Mars but its growth was arrested early by inward then 
outward migration of Jupiter. Further constraints bearing on this issue can be obtained 
by examining the conditions necessary for giant planet formation. If an embryo reaches a 
critical core mass, it cannot sustain an static atmosphere an d starts rapid gas accretion 
to fo rm a gas giant planet, such as Jupiter and Saturn (e.g 



Mizuno 



1980 



Ikomaetal 



20001 ). However, before the embryo attains the critical core mass, coUisional fragmentati on 



induced by g ravita tional stirring stalls the growth of the embryo (JKobayashi et al. 



Inaba et al. 



2010|). 



(120031 ) suggested that fr agment accreti o n enh anced by tenuous atmospheres 



allowed formation of massive cores. 



Levison et al. 



( I2OIOI ) meanwhile pointed out that 



combination of gas drag and stirring by an embryo inhibits the accretion of fragments with 
radii larger than 30 m onto the embryo. However, the typical fragment radius is smaller 
than 30 m and hence t he obstacle to fragment accretion is negligible (see the discussion in 



Kobayashi et al. 



20101 ). Since the removal of fragments due to radial drift is more effective 
than accretion, embryos cannot reach the critical core mass starting from small initial 
planetesimals. Although the core formation must be completed before gas depletion, the 
formation timescale of a core is longer than r^ for large planetesimals, for which initial 
runaway growth or cores is slow. Therefore, Jupiter's core is likely to have formed from 
accretion of intermediate size planetesimals in a massive disk (.Kobayashi et al.ii201ll ). 



Jupiter's gas accretion was truncated by its gap opening. At the edge of the gap, a radial 
pressure maximum was created. Small fragments accumulated at the pressure maximum , 



which might have lead to the formation of Saturn's core in ~ lO^yr (JKobayashi et al. 



2OI2I ). Heat released by solid accretion stabilizes atmospheres of cores and delays the 
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rapid gas accretion to for m gas giant planet s, and hence high sohd accretion rates of cores 



result in large cores (e.g., 



Ikoma et al 



Jupiter's core (ISaumon and Guillot 



formation models. 



20001). Saturn's core is estimated to be larger than 



20041. which is difficult to explain in classical core 



Lambrechts and Johansen 



(120121 ) recently proposed a core formation 



scenario by accretion of sub-meter sized pebble, which also produces a higher accretion 
rate for Jupiter unless the pebble surface density has a fiat radial profile (the radial slope 
(ilnSs/cilna > —0.5). A promising model to explain Saturn' s large core relative to Jupiter 



is if Saturn's core was formed at the gap opened by Jupiter (JKobayashi et al 



20121 ) 



In the context of the model of iKobayashi et al.l (120 111 . |2012| ). which uses the same 
prescriptions as used in the present contribution, Jupiter and Saturn can be formed in 
limited conditions for planetesimal sizes and disk masses. Large planetesimals tq ^ 30- 
200 km i n a massive disk a: > 3 are required to form a massive core of SM^ at 5.2 AU for 



Jupiter (IKobayashi et al. 



201ll ). On the other hand, the larger core of Saturn requires very 



rapid formation, which is realized by fragment accretion at the pressure maximum caused 
by Jupiter's gap opening in the solar nebula. T his yields an additiona l constraint described 



empirically by the relationship x > a/^o/I km ( Kobayashi et al. 



2OI2I ). Both conditions for 



formation of the cores of Jupiter and Saturn are satisfied for x > 7 and tq ^ 30-100 km (see 

Fig. ED. 

As shown in Figs. [3H6|, the conditions required to explain the accretion of the cores of 
Jupiter and Saturn are inconsistent with those inferred for Mars. The condition required 
for Jupiter and Saturn formation produces Earth-sized or larger embryos in the terrestrial 
planet region and the predicted e^^^VFmantie adopting /^nantie — 12 for E arth's mantle is 



> 10, which is much larger than e Wma.nt]e , 



Kleine et al. 



2009; 



Dauphas and Chaussidon 



1.85 in terrestrial rocks ( jJacobsen 



201 ih . 



2005; 



One way to try to reconcile the apparently inconsistent conditions inferred for Mars 
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on the one hand and Jupiter and Saturn on the other hand is to modify the disk surface 
density. In the discussion above and in Fig. ISHHl we assumed Sg oc Sg oc a^^-^. It is possible 
for instance that the disk surface density had a shallower radial slope of Sg oc Ss oc a~^. 
The conditions for forming the cores of Jupiter and Saturn become x > 2.4 and tq > 40 km. 
If we use a massive disk of 0.1 solar mass required by Jupiter formation, the radius of initial 
planetesimals for Mars formation becomes less than 5 km for Sg oc Sg oc a^^. Changing the 
disk density profile does not seem to solve the discrepancy that exists between Mars and 
the giant planets. 

A way to reconcile the constraints derived from inner and outer solar system bodies 
is if planets formed in a relatively massive solar nebula (x > 5) but planetesimals sizes 
were different in the inner and outer regions. In this scenario, Jupiter and Saturn cores 
would have formed from planetesi mals of 50 km radius w hile Mars formed from smaller 



planetesimals of < 10 km radius. 



Morbidelli et al. 



( 120091 ) proposed that 100 km or larger 



pla netesimals were necessa ry to explain the mass distribution of asteroids in the main belt, 
but IWeidenschillingI (120111 ) concluded that planetesimals of 100 m size could also reproduce 



the asteroid mass distribution. T 
mechanisms {e.g. 



Okuzumi et al. 



le sizes of planetesimals ar e related to their formation 



2012; 



Johansen et al. 



20061 ). The heliocentric gradient in 



planetesimal sizes discussed here may be related to the inferred contrast in solid surface 
densities and material properties across the snow line. 
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A. Derivation of £:^*^W evolution 

The mantle ^^^W value evolves due to extraterrestrial delivery, removal of metal to 
core, and radioactive decay of ^^^Hf. A planetary embryo mainly grows through collisions 
with small bodies. If the sizes of bodies are much smaller than the thickness of the embryo 



mantle, the sinking meta 



breaks up in tiny 



2003 



droplets tha t can equilibrate isotropically with 



Samuel 



20121 ) ■ On the other hand, one third of 



the surrounding mantle ( iRubie et al. 

embryo mass growth is prov ided by collisio ns with similar rnass bo dies ([Chambers 



Kobayashi et al. 



2010 . 



201ll ). According to 



Morishima et al 



2006 



( 120121 ). contribution of large 



bodies to the e^^^W evolution of the martian mantle is negligible. Therefore, we assume 
that the bodies accreted by Mars are fully equilibrated with the mantle before metal 
removal into the core. Accordingly, the equation governing the ^^^W isotopic evolution is. 



d 

"IT (^mantle 



ri82 



W] 



mantle J 



[^'^W]cHUR^ - D"^ rw] 



dMr, 



dt 



mantle " 



dt 



+ A [ Hi J mantle ^mantle) 

(Al) 



where Mmantie and Mcore are, respectively, the mantle and core masses (M = Mmantie + ^core), 
D"^ is the W metal/silicate partition coefficient. The mass balance for W in the bulk planet 
can be written as. 



[W]mantleMmantle + ^"^[Wj^antleM 



[W] 



CHUR 



M. 



(A2) 



On the other hand, Hf does not go to the mantle and the mass balance is given by 



[Hf] mantle ^mantle — [Hf]cHUR^- 



(A3) 



We assume that the core mass fraction, 7 = Mcore/M, is constant. We then obtain, from 
Eqs. ^M) and flA3l) . 

,Hf/W _ (Hf/W)mantle . _ ^^^7 ^ , ,. 

/mantle = TTTTTTTT^ -"- " I ~- K^^J 



mantle (Hf/W)cHUR ^-7 

The variation of R^ = C'^W/''''WU..ntie and i?cHUR = (''W/'''W^)chur is given by 
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Eq. dM]) divided by ^^^WmantieM using Eqs. ( R2l) : 

(1-7) i?m ,, + —3— = (1 - 7 + I^ 7)^CHUR- 



dt dt J dt 



"'' V ^^ / mantle 

Using Eq. (jUD and e'^m/^^'WU^^,i, = (1 + /Zntic)('''Hf/i82w)cHUR and dRcnvn/dt = 
A(i82Hf/i84w)^j^^j^^ we then have 

— (i?m - ^CHUr) = -(/mantle + l)(^m - ^CHUr)— ^^ ^ A/^^^^j^ ( ^^^ j . (A6) 

\ / CHUR 



Integration of Eq. (]A6P results in 



, , „Hf/W 



CHUR,0 ' 

Applying the definition of e^^'^W, Eq. (IA7p can be rewritten as Eq. 
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Figure Caption 

Fig. [T] — (a) Time evolution of embryo mass at 1.5 AU in the MMSN disk {x = 1), 
starting from different planetesimal radii tq. The gray dotted line indicates the mass 
of Mars, (b) Evolution of e^^^Wmantic obtained from the mass evolution shown in panel 
(a). The filled c ircle and the error bar represent the ^^^W excess in the martian mantle. 



e^^^^WMars mantle (jKleine et al. 



2009 



Dauphas and Pourmand 



201l|). 



Fig. |2] — Same as Fig. [T]but for x = 2.7. 

Fig. in] — Constraints from final embryo masses and e^^^W on the initial radius of 
planetesimals, tq, at 1.5 AU (a) and 1 AU (b) in the xxMMSN disk. Embryos reach the mass 
of Mars on the solid curve. The region between the thin curves produces embryos within a 
factor 2 of the mass of Mars. On the thick dashed curve, the predicted e^^^Wmantio values 
correspond to the mean value obtained from Martian meteorites. The thick dashed curve 
indicates the parameter space where e^^^Wmantie.up = ^^^^^Wmantie + Ae-'^^^Wmantie = 2.49. 
Below this limit, the accretion timescale of Mars would be too long to explain e^^^Wmantie 
measured in Martian meteorites. The colored region at small tq represents the condition 
satisfying the mass and e^^^Wmantie value of Mars and the other at Tq ^ 40-100 km indicates 
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the parameters for the formation of Jupiter and Saturn cores. For reference, the right 
vertical axis represents the disk mass corresponding to x in the left axis, assuming that the 
inner and outer edges of the disk are at 0.4 and 30 AU, respectively. 

Fig. m — Same as Fig. [3^, but including the accretion enhancement by planetary 
atmosphere (a) and gas drag (b). 

Fig. M — Same as Fig. [3^, but for different gas dissipation timescales t^ = 5 Myr (a) 

and IMyr(b). 

Fig. M — Same as Fig. [Sb (at 1 AU), but introducing the isolation time of mars, 
Ti = 20 Myr (a), 5 Myr (b), after which Mars does not grow at all. 
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